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The interplay between traffic dynamics and epidemic spreading on complex networks has received increasing 
attention in recent years. However, the control of traffic-driven epidemic spreading remains to be a challenging 
problem. In this Brief Report, we propose a method to suppress traffic-driven epidemic outbreak by properly 
removing some edges in a network. We find that the epidemic threshold can be enhanced by the targeted cutting 
of links among large-degree nodes or edges with the largest algorithmic betweeness. In contrast, the epidemic 
threshold will be reduced by the random edge removal. These findings are robust with respect to traffic-flow 
conditions, network structures and routing strategies. Moreover, we find that the shutdown of targeted edges 
can effectively release traffic load passing through large-degree nodes, rendering a relatively low probability of 
infection to these nodes. 
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Human society has always suffered from various viruses, 
such as AIDS, H1N1 influenza and computer virus. As 
the rapid development of complex network theory JIl 1. 
much effort has been dedicated to understand dynamical pro¬ 
cesses of epidemic spreading on complex networks in the past 
decade y(-[l 1 ]. Propagation is usually assumed to be driven by 
reaction processes, in the sense that every infected node trans¬ 
mits diseases to all its neighbors at each time step, producing 
a diffusion of the epidemics in the population. However, in 
many realistic situations, even when there is a link connect¬ 
ing two nodes, infection will not propagate unless some kind 
of traffic happens between the nodes. For example, a com¬ 
puter virus can spread over Internet via email-exchanges. In 
the absence of such data packet transmission, even if there is 
a path linking two computers, an infected computer will not 
be able to infect the other one. Another example is that air 
transport tremendously accelerates the propagation of infec¬ 
tious diseases among different countries. 

The first attempt to incorporate traffic into epidemic spread¬ 
ing is based on metapopulation model Il2l4l911 . This frame¬ 
work describes a set of spatially structured interacting sub¬ 
populations as a network, whose links denote the traveling 
path of individuals across subpopulations. Each subpopula¬ 
tion consists of a large number of individuals. Recently, Mel- 
oni et al. proposed another traffic-driven epidemic spread¬ 
ing model (20], in which each node of a network represents 
a router and the epidemic can spread between nodes by the 
transport of information packets. 

One of the most important issues in the study of epidemic 
spreading is how to control the prevalence of infection. To 
suppress the traffic-driven epidemic spreading, a variety of 
strategies have been considered, such as the restriction of 
traffic flow na, the selection of routings |2^] and heteroge¬ 
neous curing rate [(23], etc. In this Brief Report, we propose 
a method to control traffic-driven epidemic spreading based 
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on edge-removal strategies. The principle of edge-removal 
strategies is to affect the spreading dynamics of epidemics (or 
virus) by deleting some edges in the underlying network. It 
has been recognized that edge-removal strategies can greatly 
influence the dynamics of synchronization 124], evolutionary 
games ll25ll and traffic ll26h . In Ref. |[27t l. Zhang et al. found 
that, both random and targeted deletion of edges can sup¬ 
press the outbreak of reaction-based epidemic. In contrast, 
we will show that random and targeted edge-removal strate¬ 
gies play different roles in the traffic-driven epidemic spread¬ 
ing. Specifically, we have found that the random shutdown 
of edges decreases the epidemic threshold, while the targeted 
shutdown of edges increases the epidemic threshold. 

Following the work of Meloni et al. [20], we incorporate 
the traffic dynamics into the susceptible-infected-susceptible 
model |[28il of epidemic spreading as follows. In a network 
of size N, at each time step, A N new packets are gener¬ 
ated with randomly chosen sources and destinations, and each 
node i can deliver at most C, packets toward their destina¬ 
tions. Packets are forwarded according to a given routing al¬ 
gorithm. The queue length of each agent is assumed to be 
unlimited. The first-in-first-out principle applies to the queue. 
Each newly generated packet is placed at the end of the queue 
of its source node. Once a packet reaches its destination, it 
is removed from the system. Nodes can be in two discrete 
states, either susceptible or infected. After a transient time, 
the total number of delivered packets at each time will reach 
a steady value. Subsequently, an initial fraction of nodes po 
is set to be infected (we choose p (J = 0.1 in our numerical 
experiments). The infection spreads in the network through 
packet exchanges. All packets queuing in an infected node 
are infected, while all packets in a susceptible node are un¬ 
infected. A susceptible node has the probability (3 of being 
infected every time it receives an infected packet from any in¬ 
fected neighboring nodes. With probability 1 — /3, the virus in 
an infected packet will be cleaned by antivirus software in the 
susceptible node. The infected nodes are recovered at rate // 
(here, we set p, = 1). 

In the following, we carry out simulations systematically by 
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FIG. 1: (Color online) Density of infected nodes p as a function 
of the spreading rate f3 for the null (/ = 0), RS (/ = 0.3), BS 
(/ = 0.02) and DS (/ = 0.02) cases. The packet-generation rate 
A = 1. Each curve is an average of 10 2 different realizations. 


employing traffic-driven epidemic spreading on the Barabasi- 
Albert (BA) scale-free networks (29] with the shortest-path 
routing algorithm |f30l [3ltl . Initially, the size of BA network 
is set to be A” = 5000 and the average degree of the network 

(k) = 10. Moreover, we assume that the node delivering ca¬ 
pacity is infinite, so that traffic congestion will not occur in 
the network. 

Three edge-removal strategies are considered respectively. 

(l) The random strategy (RS): We randomly remove a frac¬ 
tion / of edges from the network. (II) The betweenness- 
based strategy (BS): We rank the edges in descending order 
according to their algorithmic betweenness. The algorithmic 
betweenness of an edge is the average number of packets pass¬ 
ing through that edge at each time step in the steady state. We 
close a proportion of edges at the top of the ranking list. (Ill) 
The degree-based strategy (DS): We define the significance 
Gij of an edge by the product of the degrees of two nodes i 
and j at both sides of the edge, i.e., G ij = hi x kj. After 
computing the significance of all edges, we rank the edges in 
descending order according to their significance. A propor¬ 
tion of edges at the top of the ranking list are removed from 
the network. For all three strategies, disconnected networks 
are avoided. 

Figure|T|shows the density of infected nodes p as a function 
of the spreading rate j3 for the null (/ = 0, i.e., no edges 
are shutdown during the epidemic spreading process), RS, BS 
and DS cases. We observe that for each case, there exists an 
epidemic threshold j3 c , beyond which the density of infected 
nodes is nonzero and increases as /? is increased. For /3 < 13 c , 
the epidemic goes extinct and p = 0. 

Figure [2] shows the ratio of /3 C (/) to /3 C (0) as a function of 
the fraction of deleted edges / for the cases of RS, BS and 
DS. Here j3 c (0) is the epidemic threshold for the null case and 
/3 C (/) represents the epidemic threshold under the condition 
that a fraction / of edges in the network are deleted. From 
Fig. El we can see that /3 c (f)//3 c ( 0) > 1 and /3 c (f)//3 c (0) in- 
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FIG. 2: (Color online) The relative ratio /3c(f)//3 c (0) as a function 
of the fraction of deleted edges / for RS, BS and DS. The packet- 
generation rate A = 1 and /3 C (0) ~ 0.043. Each data point results 
from an average over 10 2 different realizations. The curves are the¬ 
oretical predictions according to Eq. (1) and Eq. (2). The solid, 
dashed, and dotted curves correspond to the theoretical predictions 
for RS, BS and DS, respectively. 


creases with the increment of / for the cases of BS and DS, in¬ 
dicating that targeted edge-removal strategies can effectively 
suppress the outbreak of epidemic. As shown in Fig. [3 com¬ 
pared with that of the null case, the epidemic threshold can be 
enhanced more than 50% when only one percent of targeted 
edges are cutting-down. It is also noted that the epidemic 
threshold in the case of DS is a little larger than that in the 
case of BS, given that the same fraction of edges are deleted. 
For RS, however, /3 c (/)//3 c (0) is found to decrease as / in¬ 
creases, demonstrating that the random edge-removal strategy 
is failed to inhibit the spreading of epidemic, but rather en¬ 
hances its propagation. 

According to the analysis of Ref. (20], the epidemic thresh¬ 
old for uncorrelated networks is 
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where b a \ g is the algorithmic betweenness of a node 1 I 32 L [330 
and (•) denotes the average of all nodes. The algorithmic be¬ 
tweenness of a node is the number of packets passing through 
that node when the packet-generation rate A = 1 /N 1132113311 . 
For the shortest-path routing protocol, the algorithmic be¬ 
tweenness is equal to the topological betweenness (6 a i g = 
6top) and (6 a i g ) = ( D)/(N — 1), where ( D) is the average 
topological distance of a network. Here, the topological be¬ 
tweenness of a node k is defined as 


Lfc _ 1 \ ' (fc) 

top N{N-l)2g <Jij 


( 2 ) 


where <Jij is the total number of shortest paths going from i to 
j, and (k) is the number of shortest paths going from i to j 
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FIG. 3: (Color online) The algorithmic betweenness b a i s (k) as a 
function of the original degree k for the null (/ = 0), RS (/ = 0.4), 
BS (/ = 0.06) and DS (/ = 0.06) cases. Each data point results 
from an average over 10 2 different realizations. 



k 


FIG. 4: (Color online) The dependence of the density of infected 
nodes pk on the original degree k for the null (/ = 0), RS (/ = 0.4), 
BS (/ = 0.06) and DS (/ = 0.06) cases. In all cases, the packet- 
generation rate A = 1 and the density of infected nodes p ~ 0.1. 
Each data point results from an average over 10 2 different realiza¬ 
tions. 


and passing through k. The average topological distance of a 
network is given by ( D ) = JA , • dij/[N(N — 1)], where dij 
is the shortest distance between i and j. Combining Eq. (1) 
and Eq. (2), we are able to calculate the theoretical values of 
0c{f)/Pc( 0). In Fig. [3 we notice that the theoretical predic¬ 
tions agree well with the numerical results. 

To show how different edge-removal strategies affect traffic 
flow on the nodes with different degrees, we display in Fig. [3 
the dependence of algorithmic betweenness b a \ g {k) on degree 
k. Here the degree of a node is calculated before the imple¬ 
mentation of deleting edges. From Fig. [3] one can see that for 



FIG. 5: (Color online) The epidemic threshold f3 c as a function of the 
packet-generation rate A for the null case (/ = 0), the RS case (/ = 
0.06), the BS case (/ = 0.06) and the DS case (/ = 0.06). For all 
cases, the delivery capacity of a node i is equal to its original degree, 
that is Ci = hi. The arrows mark the critical packet-generating rates 
A c . For the null case (/ = 0), A c « 0.7; For the RS case (/ = 0.06), 
Ac fa 0.6; For the BS and DS cases (/ = 0.06), A c fa 1.4. Each data 
point results from an average over 10 2 different realizations. 


both the null and RS cases, b a \ g {k) increases as the increas¬ 
ing of fc, and the relationship between b a \ g {k) and k follows 
a power-law form as b a \ g (k) ~ k v . The exponent v is almost 
the same for the null and RS cases. In addition, we can also 
observe that, for large values of k, b a \ g {k) is much smaller in 
the cases of BS and DS as compared to that in the cases of null 
and RS. This point is understandable, since the targeted dele¬ 
tion of edges makes many transport paths bypass large-degree 
nodes and reroute via moderate-degree nodes, hence decreas¬ 
ing the algorithmic betweenness of those hub nodes. Conse¬ 
quently, as shown in Fig. [3 the highest values of b a \ g {k ) are 
refered to those medium-degree nodes in the cases of BS and 
DS. 

We define pi~ as the density of infected nodes of degree k. 
Figure |4] features the dependence of pt on k for the null, RS, 
BS and DS cases. Combining Figs.[3and 0] we can observe 
that the algorithmic betweenness is positively correlated with 
the risk of being infected. As shown in Fig. [3 pk increases as 
k increases for the null and RS cases. Compared with these 
two cases, the probability of being infected for large-degree 
nodes is greatly reduced in the cases of BS and DS. 

We now turn our attention to a more realistic situation 
where the node delivering capacity is finite. The main dif¬ 
ference with the infinite-capacity case is the possibility of the 
emergence of traffic congestion in the network, which oc¬ 
curs when the packet-generating rate exceeds a critical value 
A c 0)39■ Specially, we set the delivery capacity of a node i 
to be equal to its original degree, that is Ci = ki. The epi¬ 
demic threshold /3 C as a function of the packet-generation rate 
A for the null, RS, BS and DS cases are depicted in Fig. [3 We 
see that /3 C decreases and stabilizes at a constant value as A 
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FIG. 6: (Color online) (a) The relative ratio /3c(/)//? c (0) as a func¬ 
tion of the fraction of deleted edges / for the RS, BS and DS cases. 
The Internet at the autonomous system level and the shortest-path 
routing are applied. The packet-generation rate A = 0.25 and 
/3c(0) ~ 0.045. (b) The relative ratio /3 c (f)/f3 c ( 0) as a function 
of the fraction of deleted edges / for the RS. BS and DS cases. The 
BA networks and the random-walk routing are used. The packet- 
generation rate A = 0.02 and /3 C { 0) — 0.058. For both (a) and (b), 
the node delivering capacity is infinite. Each data point results from 
an average over 10 2 different realizations. 


increases. We also observe that, the random-deletion strategy 
cannot increases the epidemic threshold while the targeted- 
deletion strategies (BS and DS) can effectively enhance the 
epidemic threshold, regardless of that the traffic is in the free¬ 
flow state (A < A c ) or in the congested state (A > A c ). 

Finally, we examine the performance of our proposed 


strategies in controlling epidemic spreading by considering 
our model on different network structures and with alterna¬ 
tive routing protocols. In Fig. [6ja), we present the simula¬ 
tion results for different edge-removal strategies on the Inter¬ 
net at the autonomous system level l34tl . where the network 
size N = 6474 and the average degree (k) = 3.88 before cut¬ 
ting edges. The packets are delivered following the shortest- 
path routing. In Fig. |6]T>), we carry out our studies on the 
BA networks, where packets are forwarded according to a 
random-walk routing algorithm, i.e., a packet is delivered to 
a randomly selected neighbor. For the random-walk routing, 
the algorithmic betweenness b a \ g of a node is proportional to 
its degree |[35l [36h - As shown in Fig. [6j the conclusion that 
the suppression of epidemic outbreak by the targeted edge- 
removal strategies and the promotion of epidemic outbreak by 
the random edge-removal strategy, is still valid. 

In conclusion, we have studied the impact of edge-removal 
strategies on traffic-driven epidemic spreading. The shutdown 
of links in terms of their algorithmic betweeness or of links 
connecting large degree nodes, are found to be quite efficient 
in suppressing epidemic spreading. Contrary to previous stud¬ 
ies on reaction-based epidemic 12711 . we find that the random 
shutdown of edges accelerates the outbreak of traffic-driven 
epidemic. Furthermore, compared to the deletion of edges 
with the largest algorithmic betweeness, the shutdown of links 
connecting large-degree nodes is proved to be more effective 
in enhancing the epidemic threshold. Thus, according to our 
present studies, the targeted link-closing method can be used 
to control the spreading of computer virus in the Internet. For 
example, we can temporarily close links between large-degree 
nodes at the time of virus outbreak, and recover these links 
after virus is eliminated from the system, which could be re¬ 
alized readily by special softwares. 
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